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We explore the fluorescence signals from a pair of atoms driven towards Rydberg states on a 
three-level ladder transition. The dipole-dipole interactions between Rydberg excited atoms sig¬ 
nificantly distort the dark state and electromagnetically induced transparency behavior observed 
with independent atoms and, thus, their steady state light emission. We calculate and analyze the 
temporal correlations between intensities and amplitudes of the signals emitted by the atoms and 
explain their origin in the atomic Rydberg state interactions. 
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I. INTRODUCTION 

The measurement of temporal field correlation func¬ 
tions plays a fundamental role in the demonstration of 
nonclassical properties of light and in the analysis of the 
underlying dynamics of quantum emitters [T]. Field- 
field correlations and intensity-intensity correlations, re¬ 
flecting wave and particle properties of light, respec¬ 
tively, have thus been studied for the emission by a sin¬ 
gle atom EHg and by few atoms [H], while particle and 
wave aspects have been jointly addressed for amplitude- 
intensity correlations in the emission by a single atom 
[THin]- While the correlation functions can be expressed 
as two-time averages in the Heisenberg picture and can 
be calculated by use of the quantum regression theorem 
[HIE], a simple intuition for their behavior can be ob¬ 
tained by alternatively considering the conditional dy¬ 
namics of the light emitting system. E.g., after a photon 
counting event, the atom is put in its ground state, and 
its subsequent transient dynamics shows damped popu¬ 
lation and coherence oscillations, which are naturally re¬ 
flected by the temporal modulation of the light emission 
characteristics. 

In this article, we consider the intensity and ampli¬ 
tude correlations in the fluorescence emitted by a pair of 
atoms that are both excited on three-level ladder tran¬ 
sitions under electromagnetically induced transparency 
(EIT) conditions [T^ . Such atoms will emit only very 
little light, and since an emission event on the short¬ 
lived lower transition is accompanied by a quantum jump 
of the atom into the ground state, both the field inten¬ 
sity and amplitude will show strong transient dynamics 
around the rare counting events [ini]. In the absence of 
interactions a pair of atoms will emit light in a mutually 
uncorrelated manner, but as shown in [15j . a strong inter¬ 
action between two atoms which are both in the upper 
exited state, leads to correlations between the emitted 
intensities. 
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The physical situation can be implemented with the 
use of highly excited Rydberg states as the upper state 
in the level scheme shown in Fig. 1. Such states have 
long lifetimes, and due to their large dipole moments, 
two nearby atoms experience a large interaction energy 
shift, which in turn shifts the resonance condition for the 
excitation of both atoms. This effective detuning may 
be large enough to prevent the excitation (the Rydberg 
blockade [THUTH] ). or it may merely detune the state with 
two excited atoms. In either case, the EIT properties are 
distorted and the atoms emit more light [niiiMii]- 

Rydberg atoms have received much recent attention 
due to their potential application in quantum information 
processing |53] where the Rydberg blockade mechanism 
can be used to mediate controlled quantum gates mM- 
EIT arises from the destructive interference between dif¬ 
ferent absorption and emission processes [13], benehcial 
for slowing of light [25] and quantum nondemolition inter¬ 
actions [35], while the Rydberg blockade in combination 
with EIT is attracting interest as it paves the road for 
strong cooperative optical nonlinearities ]Il|3ni|31|37]- 

130]. 

We shall characterize the field correlations from a pair 
of atoms, exploring how they may on the one hand serve 
as a probe of the Rydberg interactions, and, on the other 
hand as a potential heralded source of nonclassical radia¬ 
tion. In Sec. II, we introduce the master equation of our 
interacting atoms and we calculate and analyze two-time 
correlation functions of the emitted radiation. In Sec. 
Ill, we turn our attention to three-time correlation func¬ 
tions, addressing in particular how the field amplitude 
emitted by one atom behaves between counting events 
and quantum jumps of the two emitters. In Sec. IV, we 
conclude and summarize our results. 


II. DENSITY MATRIX AND TWO-TIME 
CORRELATIONS 

In this section we analyze the two-time correlations 
using the master equation and quantum regression the¬ 
orem. In Sec. m we present the master equation 
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FIG. 1. (Color online) Level diagram of two three-level ladder 
atoms with a lower state 11 ), a short-lived intermediate state 
12), and a long-lived Rydberg state |3). The atoms are coupled 
to each other by the dipole-dipole interaction V 12 between the 
Rydberg states. In both atoms the lower transition, with a 
line width 71 due to atomic decay, is driven by a resonant 
probe field with Rabi frequency Hi, and the upper transition 
subject to decay 72 and dephasing 7 ph of the upper level, is 
coupled to a resonant laser field with Rabi frequency Q. 2 - 


of two three-level ladder atoms subject to Rydberg- 
Rydberg interactions. The intensity-intensity correla¬ 
tions and intensity-amplitude correlations are discussed 
in Sec. m and in Sec. |IIC[ respectively. 


A. The physical system and the master equation 

Our physical system involves two atoms excited by a 
pair of laser fields towards a Rydberg state in the ladder 
configuration, shown schematically in Fig. The time 
evolution of the atomic system can be described by the re¬ 
duced atomic density matrix p, which in the Schrbdinger 
picture obeys a linear master equation. 


P = 


( 1 ) 


where 

£p= i^[iY,p] i . (2) 

j,k 

H is the total Hamiltonian of the system and takes the 
form 


H = Hx®l2+h®H2+HM, (3) 

with single-atom Hamiltonians Hj {j = 1, 2) 

H, = + n2ag^ + H.C.), (4) 

and the dipole-dipole interaction H^d between the two 
Rydberg states 


(5) 


Here = \k)jj{l\ are the atomic operators (j = 1,2; 
k,l = 1 — 3), of the jth atom, are the laser Rabi 
frequencies, and V 12 is the strength of the dipole-dipole 
interaction. The quantum jump operators Cf. of the 
jth atom account for dissipative couplings to the envi¬ 
ronment of the system, and take the form 
(i) (i) 

and C 2 = '''"ith the respective decay rates 71 ^ 2 , 

and ~'^22 a dephasing rate 

7 ph of the high-lying Rydberg state with respect to the 
two lower states. 

We recall that in the absence of decay and dephasing of 
the upper level |3), a resonantly driven three-level ladder 
atom has a dark eigenstate, |D) = j^(ng|3) — H 2 II)), 

where Hr = 7 /IH 1 P -|- IH 2 P. The system evolves into 
this state, which has no component of the intermedi¬ 
ate, short-lived state, and hence emits no photons. For 
two such atoms close enough to each other, the Rydberg 
blockade provides a detuning of the transitions that de¬ 
teriorates the dark state of the atoms. Rather than pop¬ 
ulating a product state of the dark atomic states, |DD), 
an Atom 1 Rydberg state component restricts Atom 2 to 
explore the resonant two-level transition | 1 ) e-)- | 2 ) and 
emit fluorescence photons, and vice versa m- Fluores¬ 
cence from a single atom is known to show antibunching, 
but the signal from the atom pair, instead, shows signifi¬ 
cant bunching m- This is qualitatively explained by the 
system exploring the eigenstates of the full Hamiltonian 
H (3), which have only single occupancy of the Ryd¬ 
berg state |3), but which have also acquired population 
of the product state component |22) [3T]. When a photon 
is emitted by one of the atoms, it heralds an increased 
probability for the other atom to be in the short-lived 
excited state and emit a subsequent photon. 


B. Intensity—intensity correlation functions 


The correlations in photon counting from atoms i and 
j are evaluated as normal ordered correlation functions of 
the atomic dipole lowering and raising operators [Din], 
and are conveniently normalized by the product of the 
mean emission rates to yield the correlation function 



{<^2lit)^22it + 'rWl2it)} 

{^kit)){<^22i't + T)) 


( 6 ) 


with i,j = 1,2. This correlation function is readily calcu¬ 
lated by the master equation and the quantum regression 
theorem dlH], and it is illustrated for a characteristic 
example in Fig. see also m- We observe that the 
cross correlation is subject to damped Rabi oscillations 
with a period of 27r/HR. This is because the quantum 
regression theorem provides a linear set of equations for 
all two-time correlators with exactly the same coefficients 
[mill] as the master equation ( 1 ), but we can also in¬ 
terpret the correlation function in ( 6 ) as a consequence 
of the measurement back action. Conditioned on the de¬ 
tection event the joint state of the two atoms undergoes 
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FIG. 2 . (Color online) Intensity-intensity cross correlation 
on the lower transitions of two interacting Rydberg atoms. 
The parameters are: fli = O.271, ^2 = Syi, V12 = ^71, 72 = 
1 X 10 “'* 7 i, and 7ph = 1 x 10 ~'‘ 7 i. 


a quantum jump by the operator a\2 at time t, and the 
subsequent t + t dependent emission rate from atom j is 
given by its transient excited state population. The for¬ 
mal equivalence of Glauber’s photodetection theory and 
the quantum theory of measurements is further elabo¬ 
rated in Ref. [14]. 


C. Intensity—amplitude correlation functions 

In [7HTD]. it was demonstrated that combined photon 
counting and homodyne detection of the light emitted 
by an atom, also shows mutual temporal correlations. 
Phase sensitive homodyne detection may have near unit 
efficiency, it holds the potential to reveal more informa¬ 
tion than photon counting about the emitter dynamics, 
and, through tomography it may be used to fully recon¬ 
struct general quantum states of light [32]. Ensembles 
of atoms with Rydberg interactions exhibit interesting 
quantum dynamics and may yield sources of nonclassi- 
cal radiation, and in this work, we will thus address the 
influence of the Rydberg interactions on the intensity- 
amplitude autocorrelation and cross correlations. 

The detection of the field quadrature variable with 
phase 9j from atom j at time t -I- t is correlated with 
the counting of a photon from atom i at the earlier time 
t as quantified by the steady state {t — )■ 00 ) correlation 
function (r > 0): 

( 1 . 5 ). Q , Re -b r) 

^ (cr*2W)R-e[(o-2i(i + T))e‘®^] 

(7) 

We shall also determine the correlations for negative 
T, z.e., for field amplitude measurements preceding the 
count events, and due to the normal ordering requirement 
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FIG. 3 . (Color online) Amplitude-intensity autocorrelation 
5 ii( 22)('''’'5'/2) (S') snd cross correlation ^^2(21) (t n/S) (b) on 
the lower transitions of two three-level atoms with excited 
state (Rydberg) interaction V12 = 0 (black solid line), V12 = 
h'fil 2 (red dotted line) and fiyi (blue dashed line). The other 
parameters are the same as in Fig. 


[Tj, the correlation function is here given by a different 
expression (r < 0): 



lim 

>-oo 


Re[(cr^i(t) e'^^cr^ 2 (f-T))] 

(CT22(^-^))R-e[(cr^i(t))eiG] 


( 8 ) 


In Fig. I^a), we show the intensity-amplitude corre¬ 
lation function in Eqs. (7) and (8) for one atom, i.e., 
i = j, in the presence of the other atom. The correla¬ 
tion function is asymmetric with faster oscillations for 
negative than for positive times r. The interpretation of 
this difference was developed for the emission by a single 
atom in [T4|: After a photon detection event, the atom 
recommences its Rabi oscillation from the ground state 
11). This is a linear combination of the dark state |D) and 
two “bright” eigenstates states separated in energy from 
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|D) by ±/iriR/2, and hence physical expectation values 
oscillate naturally at the frequency riR/2. By contrast, 
to be consistent with the emission of a photon at time t, 
the evolution prior to the photon count event correlates 
with the excited state |2) which has a vanishing overlap 
with the dark state and is thus composed only of the 
bright states separated by This causes the atomic 
observables and hence the field amplitude to oscillate at 
the higher frequency Or for r < 0. See El for a more 
precise formulation of the above retrodictive argument in 
terms of the past quantum state formalism |33j . 

In our calculations, the second atom perturbs the first 
atom by the Rydberg interaction, and we see that the 
amplitude of the oscillations of the correlation function 
is reduced when V 12 is increased, mainly because of the 
normalization by the increased mean intensity. Part (b) 
of Fig. 3, shows the behavior of the cross correlation 
function between photon counting of the emission by one 
atom and homodyne detection of the emission by the 
other atom. In the absence of interactions, there should 
be no such correlation, and the function plotted should be 
identical to unity for all times, as shown by the black solid 
curve in Fig. 3(b). This is clearly not the case for finite 
V 12 , and from our analysis of the antibunching in Fig. 2 
we understand, how detection of a photon from the first 
atom, leaves the second atom with a finite population in 
the state |2), which gives rise to the oscillatory evolution 
of its observables at frequency Hr— for both positive and 
negative r. The signal amplitude shows a modulation at 
nR/2, which we ascribe to a finite dark state amplitude of 
the atom and due to the interaction with the transiently 
evolving Rydberg state population of Atom 1. 


III. THREE-TIME CORRELATION 
FUNCTIONS 

In the previous section, the correlations between a field 
amplitude signal and earlier and later counting signals 
was qualitatively explained by the conditioned evolution 
of a monitored quantum system. In this section we elab¬ 
orate on the formalism of conditional dynamics of a light- 
emitting system and we review a recently developed the¬ 
ory of past quantum states to interpret our numerical 
results for three-time intensity correlations and intensity- 
amplitude-intensity correlations. 


A. Conditioned dynamics and the past quantum 
state 

When a system is observed continuously or at selected 
instants of time, the density matrix at time t, Pc{t) is con¬ 
ditioned on the outcome of all measurements performed 
prior to t. The solution to Eq. (1) which describes the 
un-conditioned evolution of a system is therefore modi¬ 
fied with the inclusion of stochastic terms to take into 
account the back action due to the random measurement 


outcomes. The immediate consequence is that the prob¬ 
ability for a given measurement outcome is conditioned 
on the earlier measurement outcomes, and it is this back 
action dynamics that manifests itself in the temporal cor¬ 
relation functions in optical detection. 

The same correlation also formally makes the probabil¬ 
ities for the early outcomes depend on the later events, 
and if a system has been monitored on a time interval 
\ti,t 2 ], any measurement at t G [^ 1 ,^ 2 ] will have out¬ 
comes with probabilities that depend on both the prior 
and posterior measurement results. The prior informa¬ 
tion is accounted for by Pc(t), while the posterior infor¬ 
mation is incorporated in an auxiliary matrix, E{t) [33j . 
which equals the identity, I, at time t2, and which solves 
an equation backward in time, which is the adjoint of the 
equation for pc- For our purpose, it is worth noticing that 
the adjoint of a quantum jump, p^ —)■ crj 2 PcO’ 2 i: i® up¬ 
ward transition transforming E to the excited state just 
prior to the detection event, E —> (j\iEa\ 2 ^ while in inter¬ 
vals with no probing (corresponding to the time between 
t and t -|- T in our correlation function calculations), E 
solves an equation like the master equation (1), but back¬ 
wards in time and with the first of the damping terms in 
(2) replaced by the Hermitian adjoint EC^\ 

A general measurement is described by the theory of 
positive operator valued measures, (POVMs), i.e., a set 
of operators {Um}, that fulfils The dif¬ 

ferent operators are associated with different outcomes, 
enumerated by the continuous or discrete index to, and 
for a given density matrix p(t), they yield the outcome 
probabilities, P{m) = Tr(Ump(t)IIm)- With the incor¬ 
poration of both prior and posterior knowledge, this ex¬ 
pression is replaced by [33], 

Pp{m)= Tr(U^pc(t)flLE(t)) ^ 

Eq. yields predictions that have been succesfully com¬ 
pared with the outcome of measurements on a microwave 
cavity field |34| and on a superconducting qubit [35]. 

B. Three-time intensity correlation functions 

When applied to intensity and amplitude measure¬ 
ments, Eq. (|^ reproduces the results that we obtained 
with the quantum regression theorem and illustrated in 
Figs. (2) and (3), but more importantly, it offers a for¬ 
mal theory from which we can also extract qualitative 
results and interpretations. We will in this section pro¬ 
ceed with an investigation of the three-time correlation 
function between a prior and a posterior counting event 
(at times t and t + T) and an intermediate intensity or 
amplitude measurement at time t + r. We were originally 
motivated to carry out this study by the apparent con¬ 
flict between the difference in the evolution frequencies 
in Fig. 3(a), and the results directed us to a number of 
interesting observations that we shall present below. 
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FIG. 4 . (Color online) <?[i2{22i) (c on the lower transitions 
of two interacting Rydberg atoms for different values of T = 
57)”^ (a), 107)”^ (b), and 157 )”^ (c). The other parameters 
are the same as those in Fig. 



FIG. 5 . (Color online) 5 i 22 { 2 ii) on the lower transitions 
of two interacting Rydberg atoms for different values of T = 
57)”^ (a), 107)”^ (b), and 157 )”^ (c). The other parameters 
are the same as those in Fig. 


Let us first consider the three-time intensity correlation 
function, 

^ (o'2ift)o'2i(^ + + T)a{^{t + T)(7\^{t)) 

(o'22W)(c^L(^ + 'r)) (c^f2(^+ ^)) 

with T (T) being the time interval between the first and 
the second (third) detector clicks, and i,j, k = 1,2. 

Figure 1^ shows 5 ( 12 ( 221 ) ^ function of r for 

different values of T. Near r = 0, the function reflects 
the antibunching of the signal from a single atom, while 
the ability of two different atoms to emit at the same 

/o'! 

time is witnessed by the finite value of 5112 ( 221 )^)- 

In comparison, 5(22(211) ^)> shown in Fig. reflects 

bunching for r = 0 and antibunching of light coming from 
the same atom at r = T. Note that the functions shown 
in Figs. 4 and 5 are not merely time reversed copies of 
each other. The oscillatory correlations in Fig. 4 are 
damped with time, while in Fig. 5, they evolve for long 
times with only little damping. 


C. Three-time intensity—amplitude—intensity 
correlation functions 

To examine further the evolution of the pair of atoms 
between the two photon counting events, we supplement 
the analysis of the three-time intensity correlations with 
a calculation of the expected outcome of homodyne de¬ 
tection of the field-amplitude between counting events at 
times t and t -I- T, characterized by 



^ Re [(o'li(f)cr^i(t -b t) e'^^crf2(f -b T)a{^{t))] 
(c^22W)(o'22(i + r))Re[(CT2i(f+ T))e'«^] 
with i,j, k = 1,2. 

We focus on the same cases of prior and posterior 
counting events from different atoms while probing now 
the field amplitude from the atom that caused either 
the first or the last count. We find that the correla¬ 
tions between the initial intensity measurement on one 
atom and the subsequent amplitude measurement on the 
same atom decay after a few oscillation periods in Fig. 6, 
while the correlations involving the amplitude measure¬ 
ment and the last counting measurement on the same 
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FIG. 6 . (Color online) on the lower transi¬ 

tions of two interacting Rydberg atoms for different values of 
T = 57)”^ (a), 107)”^ (b), and 207)”^ (c). The other parame¬ 
ters are the same as those in Fig. 



Y^T 

FIG. 7 . (Color online) <7^22^(211) (g^)'^/ 2 ) on the lower transi¬ 
tions of two interacting Rydberg atoms for different values of 
T = 57)”^ (a), 107)”^ (b), and 207)”^ (c). The other parame¬ 
ters are the same as those in Fig. 


atom are undamped as a function of r in Fig. 7. 

The persistent oscillations in /2) over 

time intervals much longer than the intermediate atomic 
state lifetime I/ 71 , indeed have a natural explanation in 
terms of the past quantum state formalism in Eq. (|^. 
The time evolution of pc between the detector click at 
t and the amplitude detection at time t -I- r is governed 
by the usual master equation, with the initial state given 
by the application of the jump (down) operator on the 
steady state density matrix. Likewise the matrix E is 
given by the adjoint master equation solution, propa¬ 
gated backwards in time from the value, post-selected by 
the counting event at time T. Just prior to this counting 
event, E is given by a jump (up) operator applied to the 
identity matrix |33j . 

Despite Pc{t) and E{t) having quite different dynam¬ 
ics, the linear set of equations for their matrix elements 
have identical spectra of complex eigenvalues Xj and their 
time evolution can be expanded on eigensolutions which 
evolve with complex exponential factors exp(Am'r) and 
exp(A„(r —r)), respectively. The predictions by Eq. ^ 
will hence be a combination of products of such factors, 
yielding r-independent terms (when m = n), and un¬ 
damped, oscillatory r-dependent terms from eigenvalue 



FIG. 8. (Golor online) The maximum (red dotted line), min¬ 
imum (green dash-dotted line), and mean (blue solid line) of 
Ali22(2ii)('r, T, 7r/2) around r = r /2 as a function of T. The 
other parameters are the same as those in Fig. 
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pairs with A„ = + i?/„ = . All complex eigenvalues 

occur in such pairs since the density matrix equations 
can be represented by real sets of equations. 

The sum of constant and oscillatory terms causes the 
oscillations in Fig. 7 to occur around constant levels that 
depend on T. We have illustrated this dependence in 
Fig. 8, where the lower (green dash-dotted) and upper 
(red dotted) curves show the range of variation of the 
field amplitude AIi22(21i)(tj 7’)7J'/2) around r = T/2 as 
a function of T, normalized by the intensity-intensity 
correlation between t and t + T, 


AMr,T,e,) = gbfc 


( 12 ) 


We observe that even for very long intervals T between 
the clicks, for one value of T, the amplitude oscillates be¬ 
tween zero and a large positive value, while for an only 
slightly larger T, the oscillations occur between a corre¬ 
spondingly large negative value and zero. These abrupt 
changes follow the periodicity of the intensity-intensity 
correlation function, shown in Fig. 2, which implies that 
they are rare and difficult to observe. 

While the persistent oscillations in Figs. 5 and 7 follow 
from the above eigenvalue argument, we need to explain 
why the oscillations in Figs. 4 and 6 are damped. The 
results in Figs. 4-7 are all obtained by solution of the 
same differential equations, but the conditioning count¬ 
ing events impose boundary conditions that strongly in¬ 
fluence the observed dynamics. After counting a photon 
from Atom 1, the density matrix factorizes into a prod¬ 
uct state p —>■ <J^i 2 cc |l)i(l| ® Pc, 2 , where the condi¬ 
tioned state pc .2 of Atom 2 depends on the steady state of 
the joint system, attained prior to the jump. Just before 
the counting of a photon from Atom 2 at time t -I- T, the 
matrix E acquires the form E —>■ cr 2 i’^^^i 2 ^ oc Ii 0 |2)2(2|, 
i.e., it factorizes into a product of the identity matrix Ii 
acting on the Atom 1 Hilbert space and the excited state 
projector of Atom 2. 

To understand why these limiting conditions lead to 
the observed correlations, we consider for simplicity the 
evolution between t and t + T m the absence of the Ry¬ 
dberg interaction. The matrices pertaining to the sep¬ 
arate atoms then evolve independently, and while the 
three terms, |l)i(l|, Pc ,2 and | 2 ) 2 ( 2 | all develop damped 
oscillatory dynamics, the identity matrix Ii is invariant 
under the backward time evolution of the matrix E. This, 
in turn, implies that according to Eq. the outcome 
probabilities for measurements on Atom 1 between 0 and 
T are given by the usual density matrix expression and 
hence they undergo the usual damping towards steady 
state. In our argument about eigenvalues, the terms in 
the evolution of p, damped as exp(AmT), are not multi¬ 


plied by the matching functions exp(A„(T — t)) because 
Ii is the zero eigenvalue solution of the backward evo¬ 
lution equation. The interaction between the atoms of 
course modihes the argument, but the solutions change 
and the factorization becomes invalid only gradually as 
the equations are integrated over time. 


IV. DISCUSSION 

In summary, we have studied photon correlations in 
signals emitted by a pair of interacting atoms. The inter¬ 
action signihcantly alters the dark states of the individual 
atoms, and the optical response is correspondingly large. 
This may pave the way for the use of photon correla¬ 
tion measurements as a probe of the Rydberg interaction 
strength, and, ultimately we may use the detection of 
fluorescence from one atom to herald the emission of a 
nonclassical state of light by the other atom. 

The correlations reveal a host of interesting results. In 
particular, we find strong correlations between counting 
events on one atom and field amplitude measurements on 
the same and on the other atom. We observe a strong 
asymmetry between negative and positive time correla¬ 
tions and ascribe this to the different boundary condi¬ 
tions for the equations of motion; An earlier photon 
count is accompanied by a jump into the ground state, 
while a later count event specifies the matrix E to pop¬ 
ulate the short-lived excited state, and these states have 
different expansions on the eigenstates of the ladder sys¬ 
tem Hamiltonian. 

Three-time correlation functions show oscillatory be¬ 
haviour as function of the intermediate time argument 
t T between counting events at times t and t + T, and 
these turn out to be undamped for times much longer 
than the lifetimes of the atomic states. This result is a 
consequence of the post selection by the last detection 
event, and it is explained quantitatively by the contri¬ 
bution of exponentially damped terms with both time 
arguments r and T — t. Interestingly, due to the condi¬ 
tions imposed by the detection events, this phenomenon 
occurs for some but not all combinations of correlation 
functions of the emission signals from two atoms. 

Our work focused on pairs of atoms, but we imagine 
that coherent ensemble emission from Rydberg blocked 
ensembles and from atoms with more complex level struc¬ 
ture (e.g., excited from a larger Zeeman and hyperhne 
ground state manifold) | Jtij . will show similar behavior 
and may lead to applications of the results presented. 
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